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Summary: We propose a novel method for the description of spatial patterns formed 
by a coverage of point sets representing galaxy samples. This method is based on 
a complete family of morphological measures known as Minkowski functionals, which 
includes the topological Euler characteristic and geometric descriptors to specify the 
content, shape and connectivity of spatial sets. 

The method is numerically robust even for small samples, independent of statistical 
assumptions, and yields global as well as local morphological information. We illustrate 
the method by applying it to a Poisson process, a 'double-Poisson' process, and to the 
Abell catalogue of galaxy clusters. 
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1. Introduction 



Statistical measures provide important tools for the comparison of inhomogeneous cosmo- 
logical models with observations of large-scale structure in the Universe. The discrimi- 
native power of this comparison depends chiefly on the statistical measure used. In the 
cosmological community the most frequently used measure was and still is the two-point 
correlation function of a point process such as the distribution of galaxies (e.g., Peebles 
1980 and references therein, Davis & Peebles 1983), or clusters of galaxies (e.g., Bahcall 
& Soneira 1983). This measure appears to be crude: if we consider a cosmological density 
field, then the two-point correlation function only embodies information contained in the 
power spectrum of this field; it does not provide any information on the distribution of the 
phases of the Fourier components of the field, and is limited to spatial regions which are 
well below currently observed large-scale structures such as, e.g., the "Great Wall" (Geller 
& Huchra 1989). As a consequence, completely different spatial patterns could display 
the same two-point correlation function (Martinez et al. 1990), i.e., no direct conclusions 
about the morphology of large-scale structure can be drawn. However, in general, the char- 
acteristics of the two-point correlation function as a small-scale measure are affected by 
the possible existence of large-scale inhomogeneities in an uncontrollable way. Moreover, 
they could depend strongly on 'biasing' assumptions about the relation between dark and 
luminous matter (Martinez et al. 1993, Buchert & Martinez 1993). 

Efforts to go beyond the two-point correlation function comprise the consideration of 
higher-order correlations (Peebles 1980 and ref. therein), integrated unnormalized pair- 
counts (Mo et al. 1992), void probability functions (White 1979), percolation analysis 
(Shandarin 1983), minimal spanning trees (Barrow et al. 1985), Voronoi foam statistics 
(Icke & van de Weygaert 1987, van de Weygaert 1991), cell-count variance (Efstathiou 
et al. 1990), topological/geometrical structure discriminators (Mo & Buchert 1990), and 
multi-fractal measures (Martinez 1991). However, statistical measures which are sensitive 
to the morphology of structures have been investigated mostly in other fields such as image 
analysis and pattern recognition (Rosenfeld & Kak 1976 and ref. therein, Serra 1982). 

Of those approaches which focus on global aspects of the matter distribution, the 
use of topological descriptors appears to be among the most effective with respect to the 
comparison of cosmological theories and observations. Melott and collaborators (see the 
review by Melott 1990 and ref. therein) have pursued this possibility and have investigated 
methods to probe the topology of large-scale structure by employing the Euler character- 
istic or the genus of iso-density level surfaces of a density field obtained by smoothing the 
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observed galaxy distribution. Recently, this method was modified by Brandenberger et al. 
(1993) into a discrete genus statistic. 

This approach translates the phenomenology of the large-scale distribution of matter 
we observe in the Universe into a quantitative statement of connectedness of these struc- 
tures: for different iso-density levels different topologies are realized by the same matter 
distribution. For high density levels isolated islands remain showing a bubble or cluster 
topology, for low density levels a connected structure occurs which has the topology of 
a honeycomb or swisscheese, intermediate topologies are named spongy, i.e., holes and is- 
lands are interlocking (compare Melott 1988, 1990). Also, the change of the topology of 
large-scale structure can be an evolutionary effect: while the early distribution of QSO's 
exhibits a "bubble" topology, the evolution of overdense structures into interconnecting 
sheets of superclusters displays the opposite topology as is seen in the evolution of pancake 
models typical for a EDM (Hot-Dark-Matter) cosmogony. 

The full morphological specification of spatial patterns requires topological as well as 
geometrical descriptors to characterize not only the connectivity but also the content and 
shape of figures. The aim of our paper is to point out that integral geometry supplies a 
suitable family of such descriptors, known as Minkowski junctionals. In a <i-dimensional 
ambient space there exist d + 1 of these functionals. We are primarily interested in the 
case d = 3, where the Minkowski functionals are related with familiar measures: covered 
volume, surface area, integral mean curvature, and Euler characteristic. 

In order to apply the functionals to point-sets representing galaxy samples, we proceed 
by decorating each point with a spherical ball. The scale-dependent morphological features 
of this coverage may then be explored by varying the radius of the balls. 

The Minkowski functionals are global and additive measures. Additivity allows us to 
calculate these measures effectively by summing up their local contributions, as in the case 
of the Euler characteristic. Moreover, the functionals are numerically robust against short- 
scale spatial irregularities in the coverage, which may arise via the variation of the ball 
radius. Finally, a further advantage is that the computation does not rely on statistical 
assumptions concerning the galaxy distribution. 

In the next section we briefly review the mathematical background. In Section 3 the 
results of sample calculations are described. Section 4 gives a summary. 
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2. Minkowski functionals 



Consider a set of points {xi/i = 1, . . . , N} which represent the positions of a sample of 
galaxies in three-dimensional Euclidean space IE 3 . In order to study the global features 
of this point set, we introduce a collection of neighborhoods by covering each point with 
a spherical ball B r (xi) = {x G E 3 /||x — Xi|| < r}. For any chosen value of the radius 
r the union of these balls specifies an ensemble of clusters, where two balls are said to 
belong to the same cluster if they are connected by a chain of intersecting balls. By 
employing the Minkowski functionals to measure content, shape and connectivity of the 
covering £>(r) = uf =1 B r (xi) under variation of r, we arrive at a quantitative operational 
definition for the scale-dependent morphological characteristics of the underlying spatial 
point process. 

Since the methods of integral geometry are perhaps not widely known among physi- 
cists, we compile in this section some pertinent facts (Hadwiger 1957, Santalo 1976, Weil 



The Minkowski functionals derive from the theory of convex sets. They generalize 
curvature integrals over smooth surfaces to the case of surfaces with singular edges and 
corners. In the example of the covering £>(r) such irregularities arise from the intersections 
of the decorating balls. 

Let K be a compact convex body in IE 3 , with regular boundary dK e C 2 and principal 
radii of curvature R\ and R 2 . The surface area F of K is given by 



since the area element df on dK is related with the area element da of the spherical image 
of dK under the Gaussian map by df = R\R 2 da. 

Consider the convex body K e parallel to if at a distance e, i.e.: 



1983). 




(1) 



K e = U xeK B e (x) . 



(2) 



Thus, the parallel body K e consists of all points within a distance < e from K. 
The radii of curvature of dK e are R± + e and R 2 + e. Therefore 



F £ = (Ri+e)(R 2 +e)da = F + 2-H-e + G-e 2 , 



Js 2 



(3) 
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where 

is the integral mean curvature and 



G 



[ — ^— d/ (5) 



is the integral Gaussian curvature of <9X. 

If V denotes the volume of K and V e that of K e , we have 



47T 

K = y+ / F e ,de' = V + F-e + H-e 2 + —e s , (6) 







which is known as Steiner's formula. 

More generally, Steiner's formula for convex bodies in d-dimensional Euclidean space 
may be written as 



Ve{K) = j^( d \w a {K)e« , 

a=0 ^ ' 



(7) 



whereby the Minkowski functionals W a (K), a = 0, . . . , d are defined. For d = 3 we find 
by comparison 

W Q {K) = V{K), 3W 1 (K) = F(K), 3W 2 (K)=H(K) , 3W 3 (K) = G(K) = Att X (K) , 

(8) 

with the constant x(K) = 1 denoting the Euler characteristic of a convex body. 

The expression (7) for V e (K) allows us to extend the definition of the functional W a 
to any convex body without requiring regularity of dK; compare Fig.l. For example, let 
us take a cube Q in IE 3 with edge length a. Its parallel volume V e (Q) may be found by 
inspection, 

VJQ) =a 3 + 6a 2 e + 12a-e 2 + 8— e 3 . (9) 

4 24 

The e-terms are the contributions from the rectangular prisms on the faces of Q, from the 
cylindrical sectors on its edges and from the spherical sectors located on its corners. 

Furthermore, Steiner's formula also yields the Minkowski functionals for lower dimen- 
sional ("improper") convex bodies embedded in E d , such as points, straight line segments, 
discs etc. . 

For example, in the case of a disk D r with radius r in IE 3 one finds 

V e (D r ) =2nr 2 e + n 2 re 2 + — e 3 , (10) 

o 
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and thus 

W (D r ) = V(D r ) = , 3Wi (D r ) = F(D r ) = 2nr 2 , 

3W 2 (D r ) = H(D r ) = n 2 r , 3W 3 (D r ) = G{D r ) = 4n X (D r ) = 4n . (11) 

We note that the two-dimensional disk enters here as a cylindrical prism with vanishing 
height; therefore its "content" F(D r ) = 2irr 2 is twice the area of a two-dimensional ball 
with radius r. If r — > 0, the disk degenerates into a point P = D r= o with the Euler 
characteristic x(-P) = 1 and W a (P) = 0, a = . . . d — 1. 

Let us now list some common general properties of the functionals W a : K, — > H. 
Here, /C denotes the class of closed bounded convex subsets of E d . 

Additivity : If K e K, is dissected by a planar cut into two convex parts such that K = 
K x Uif 2 , K X ,K 2 E /C, then 

WaiKi U JT 2 ) = W a (K t ) + W a {K 2 ) - W a (K! n K 2 ) , (12) 

holds for each a = 0, . . . , d. 

This relation is easily verified if K = B r (x) is a ball in IE 3 which is cut into two hemispheres 
K\ and K 2 . For instance, the surface area (4nr 2 ) of K\ U K 2 is obtained by adding up the 
surface areas of the hemispheres (2 • (2irr 2 + rrr 2 )) and subtracting the area (2yrr 2 ) of the 
disk D r = K x n K 2 . 

Motion invariance : Let Q be the group of motions (translations and rotations) in E d . The 
transitive action of g G Q on K G /C is denoted by gK. Then 

W a (gK) = W a (K) ;a = 0...d , (13) 

i.e., the Minkowski functionals of a body are independent of its location in space. 

Continuity : If K n — > K for n — > oo, K n , K E JC (with convergence defined in terms of the 
Hausdorff metric for sets), then 

W a (K n )^W a (K) ;a = 0...d . (14) 

Intuitively, this continuity property expresses the fact that an approximation of a convex 
body by convex polyhedra K n , for example, also yields an approximation of W a (K) by 
W a (K n ). 

So far we only considered the Minkowski functionals restricted to the class /C of convex 
bodies. Since we want to apply these measures to investigate the spatial pattern of cover- 
age models as described at the beginning of this section, the definition of the Minkowski 
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functionals must be extended to non-convex sets such as B(r). This extension may be 
achieved in the following way (Hadwiger 1957). 

Let 1Z denote the class of subsets of E d which can be represented as a finite union of 
sets from /C, i.e., A e TZ if and only if A = \jf =1 Ki, N < oo, Ki G /C. The class 1Z also 
includes the empty set 0. In the first step, the Euler characteristic is introduced by 

*w-{J : A 6 (15) 

and extended to 1Z via additivity, 

X (A UB) = X (A) + X (B) - X (A n 5) (16) 
for any A.BeIZ. In particular, 

xCuj^^i) = J2 x( K *) - E n + . . . + (-i)" +1 x(tfi n k 2 n . . . n k n ) , (17) 

i i<j 

which follows from eq. (12) by induction. The right hand side of eq. (17) only involves 
convex sets and may be applied together with eq. (15) to compute x(A) for any A e 1Z, 
as illustrated in Fig. 2. We also note that x : 7?. ^ ^ is motion invariant and it can 
be shown to agree with the Euler characteristic as defined in algebraic topology. In the 
cosmology literature (e.g. Melott 1990) the Euler characteristic is usually associated with 
the surface dA of a body A. For closed (d — 1) -dimensional surfaces in E d one has the 
simple relationship 

X(0A) = X (A) (1 + (-If- 1 ) . (18) 

As a second step of the extension, the Minkowski functionals are defined for A e TZ by 
W a (A) = J X (AnE a )dfi(E a ) , a = 0,...,d-l, 

(19) 

W d (A) = u dX (A) , u d = n d / 2 /T(l + d/2) . 

Here, E a is a a— dimensional plane in lE d . The integral runs over all positions (induced by 
translations and rotations) of E a , weighted with the so-called kinematical density d^{E a ) 
(Hadwiger 1957, Santalo 1976) which is related to the invariant Haar measure on the 
group of motions Q and is normalized such that for a rf-dimensional ball B r with radius r, 
W a (B r ) = u)dr d ~ a . The volume of the unit ball (B r= i) is uj d - 

In the case of convex A e JC the eqs. (19) reproduce the Minkowski functionals 
obtained previously from V e (A) (e.g. Miles 1974). 
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According to their definition (19), the Minkowski functionals on 1Z inherit from the 
Euler characteristic the property of additivity, i.e., 

i i<j 

as well as motion invariance. These features together with their "conditional continu- 
ity" (the W a are continuous when restricted to JC) specify the Minkowski functionals as a 
distinguished family of geometrical and topological descriptors: There is a remarkable the- 
orem (Hadwiger 1957) which asserts that any additive, motion invariant and conditionally 
continuous functional T on subsets A C E d , A e 1Z, is a linear combination of the d + 1 
Minkowski functionals, 

d 

^(A) = J2c a W a (A) , (21) 

a=0 

with real coefficients c a independent of A. Important consequences of this theorem are the 
"principal kinematical formulae" , which may be written concisely in the form 

M a (A n gB)dg = (fy M a .p(B)M p (A) 

M a (A) = ^—^-W a (A) , a = 0,...,d. (22) 

The integral is performed with respect to the invariant Haar measure dg of Q and runs 
over all motions of the set B, with A, B e TZ. In the case B = B e (x) and A = K e JC the 
kinematical formula for a = d reproduces Steiner's formula (7), 

/ X (KnB £ (x))d d x = V e (K). (23) 

The kinematic formulae are useful mathematical tools in stereology and stochastic 
geometry. The Minkowski functionals can be calculated efficiently for any given cover- 
age without requiring statistical assumptions about the underlying point set. However, 
when the latter constitutes a realization of a random spatial process (which is the case in 
cosmological models for structure formation) then one faces a difficult problem: to infer 
the statistics of random geometrical patterns from a generally limited amount of point set 
data. It turns out that the mean values of the Minkowski functionals provide unbiased 
stereological estimators (Weil 1983). Moreover, these mean values can be calculated ex- 
actly for the classical Boolean model (Mecke and Wagner 1991), where the centers of balls 
are distributed in E d according to Poisson's law, which is often employed as a reference 
model. 
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For convenience, we recall the result for these mean values per ball B r and for d = 3: 

<^>=I(l- e -") , (24a) 

< L >= Aitr 2 e- V , (246) 

< f >= 47rr(l - f^)e-" , (24c) 

< ^ >= 4tt < X >= 4tt(1 - Zv + ^v 2 Y- v , (24d) 

with v = 47rr 3 p/3 and the number density p of balls. 

In the general case where the statistics of a point set in a domain C E 3 with volume 
is specified by a sequence of density correlation functions {p n (xi, . . . , x n ) \n = 1, 2, . . .}, the 
mean values for the augmented coverage per unit volume are obtained from the additivity 
relation (20) in the form 

< ^ w ^)]>=Y J —^ r - J ■■ J^Pn(yi,...y n )W ot [nUBr(y l )]d s y 1 ...d 3 y n . 

(25) 

Obviously, the Minkowski functionals embody information from every order n of the cor- 
relations. 
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3. Examples 



In this section we calculate the Minkowski functionals for three examples of point processes 
in E 3 : a Poisson process, a 'double-Poisson' process and the ACO-catalogue of galaxy 
clusters. 



3.1. Calculation method 

The Minkowski functionals W a [S(r)] for a covering £>(r) = ufL 1 B r (xi) of a given point set 
{x\, . . . , xn} C E 3 may be calculated straightforwardly via the additivity relation (20). 
However, in the cases a > 1 this algorithm becomes inefficient when the amount of overlap 
between the augmented balls is excessive, since one has to compute many redundant and 
mutually cancelling terms. Therefore, we proceed, for a > 1, alternatively as follows 
(Mecke 1993): 

The covering B(r) divides the whole space by piecewise smooth spherical contours. These 
pieces are joined along circular arcs where two balls intersect, and at singular points 
(vertices) which have equal distance to the centers of three balls. Configurations where 
four or more arcs have a point in common carry negligible weight and can be ignored. 
Consider a ball B r (xi) which is only partially covered by other balls. The area Fi and the 
integral curvatures Hi and Gi are given by 

F i= [ df = 3W l (S i ) , H i = -F i = 3W 2 (S i ) , G { = ^ = 3W 3 (Si) . (26) 

An uncovered arc Aij which length 1^ on the intersection B r {xj) H B r (xj), \\xi — Xj\\ = 
ciij < 2r, gives rise to the curvature contributions 

Ha = lotijiij = 3W 2 (A ZJ ) , dj = 7 p^L == = 3W z {A ij ) , (27) 
2 ryH - {aij/2) 2 



where denotes the angle between the normals of B r (xi) and B r (xj) along 
The Gaussian curvature Gijk = 3Ws(Pijk) at an uncovered vertex Pijk on the common 
intersection of three balls centered at Xi, Xj, Xk equals the solid angle spanned by the three 
normals of the balls at and is obtained from PHuilier's formula 

[tan(-Gijfc)] 2 = 

tan( ) tan( ) tan( ) tan( ) , 
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Finally, the total values of W a [B(r)] are given by 

W 1 [B(r)]=J2w 1 (S i ) , 

i 

W 2 [B(r)] = J2w 2 (S l )-J2W2(A lJ ) , 

i ij 

W 3 [B(r)] = W 3 (Si) ~J2 W ^j) + J2 W ^ P ^ = ^x[B(r)] , (29) 

i ij ijk 

where the sums run over all uncovered surface pieces, arcs and vertices. 

In the figures below we display the reduced values of the Minkowski functionals $ Q = 
W a [B(r)]/NW a (B r ), where N is the number of sample points, and W a (B r ) = ^L r d ~ a , d = 
3; in particular, $ 3 = x[B(r)]/N. 

As a consequence of the additivity of Minkowski functionals, the statistical errors are 
always smaller than the graphical symbols used in the figures. The computation time on 
a local PC requires a few CPU-seconds for a sample of hundred points, increasing with 
N log(iV) with the number iV of points. 

In the equations (29) the global Minkowski functionals are expressed in terms of local 
measures. The latter also allow us to introduce individual Minkowski functionals W a (xi; r) 
for each ball B r (xi)(=: Bi, for short), in the following way: 

w a {x % -r) = w a (Si) -\Y, WaiA^xiBi n B,) + ^ ^ W a {P ijk ) X {Bi n B 3 n B k ) . 

j,j^i j<k,j,k^i 

(30) 

The Euler characteristics in equation (30) equal unity if the convex intersections are non- 
empty and vanish otherwise. The factors 1/2 and 1/3 arise, since an arc and a vertex are 
shared by two and three balls, respectively. Moreover, W\(Aij) = Wi(Pijk) = ^(-Fyfc) = 
0. These individual local measures (30) provide means for defining subsamples by selecting 
for fixed r those balls whose associated measures attain values in a prescribed interval. 
Examples are shown in Fig. 6 below. 
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3.2. Poisson process 

In the first example the coverage £>(r) is based on 10 4 points, independently distributed 
in a unit cube with number density p. We use the length scale x = p 1 l' i r, normalized by 
the mean distance of the points; thus, the edge-length of the cube is x cu ^ e = 21.6. The 
radius of the decorating balls is varied within the interval < x < 1.6; its maximum value 
is small enough compared with x C ube 

to avoid finite-size effects. Fig. 3a displays the values 
of $ a (x), 1 < a < 3. The curves are the mean values < <& Q (:r) > as obtained from eqs. 
(24). For small radii, B(r) consists of isolated balls with negligible overlap; therefore, each 
measure starts out at unity. As the overlap increases, the measures decrease together with 
the reduced total surface area $i(x) of the coverage. $2(3^) changes its sign, since the 
curvature of the singular edges at intersections is negative and starts to dominate the posi- 
tive contribution from the spherical parts of the surface as the overlap increases. Negative 
values of the Euler characteristic $3(2;) are typical for a highly connected, "sponge-like" 
structure with many tunnels. We also note that the first zero of $3(2;) at x p 0.44 
provides an estimate for the percolation threshold (Mecke & Wagner 1991). 



3.3. 'Double-Poisson' process 

This process was implemented as follows: 50 centers of spherical shells were randomly 
placed in a unit cube. Then, 200 points were distributed independently within each shell. 
The resulting total number of points in the cube turned out to be 6841. The edge-length 
of the cube was £ C ube = 19.0, the inner and the outer radius of each shell has been chosen 
as Xi = 3.4 and x a = 3.8. 

Fig. 4 presents the projections of 750 points from a thin planar slice cut out from the 
cube. For comparison, we also show a two-dimensional Poisson process with 500 points. 

The measures $> a (x) for the 'double- Poisson' process are given in Fig. 3b. In com- 
parison with Fig. 3a, the x— range with negative Euler characteristic is extended: As x 
increases, the balls start to overlap earlier and they fill the space later. This feature is also 
reflected in the filaments and bigger voids recognizable in Fig. 4. The minimum value of 
the Euler characteristic is higher than in the pure Poisson case, indicating that now the 
structure is less "spongy" with fewer (but bigger) tunnels. 
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3.4. AC O- catalogue 

For the third example of a point process we take the ACO-catalogue of galaxy clusters 
(Abell et al 1989, see also Postman et al 1992 and the recent paper on the topological 
characterization of the catalogue by Rhoads et al. 1993). We consider a spherical section 
of the northern hemisphere with an opening angle of 120° and a radius of 450 h _1 Mpc 
or z = 0.15. Each cluster in this section is taken into account, but we neglect thinning 
in distance and other possible deviations from a 'fair sample'. In this way we obtain 923 
clusters of galaxies; the radius of the section is auction = 9.6 in normalized units, again big 
enough to avoid edge effects, as in the previous examples. Since the clusters are already 
extended objects, the ball radii are chosen such that x > 0.06 (r c i uste r > 3 h _1 Mpc). In 
Fig. 3c the Minkowski measures for the coverage of this observed point process are plotted. 

To exhibit the significance of these measures we display in Fig. 5 the data for the 
Euler characteristic per unit surface area x*( x ) '•= x( x )/F(x) for the three examples. The 
enhancement of x*( x ) a t larger radii results from the reduction of the size of uncovered 
tunnels. The similarity between the data obtained from the ACO-catalogue and those 
from the 'double-Poisson' process is remarkable in view of the fact that the parameters 
chosen in the latter have not been adjusted or optimized. 

In Fig. 6 we show planar projections of three-dimensional ACO-subsamples selected 
by the values of the individual Euler characteristics of partially covered balls, as described 
in Section 3.1. Only the centers of the balls are displayed. 
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4. Concluding Remarks 



We have presented a method to analyze the spatial pattern of galaxy distributions by 
employing the Minkowski functionals which provide a complete family of morphological 
measures. In order to apply this concept to the point-set representing a sample of galaxies, 
we introduced a coverage by marking each member of the point-set with the center of a 
ball. The common radius of the balls is used as a diagnostic parameter for displaying the 
morphological features of the coverage on different length-scales, i.e., with varying spatial 
resolution. In three-dimensional space the family of Minkowski measures consists of the 
Euler characteristic, the covered volume, the surface area of the coverage, and its integral 
mean curvature. 

The important aspects of the method may be summarized as follows. 

• On the practical side, the calculation of the Minkowski functionals is precise and CPU 
time effective. In contrast to the currently popular approach (see, e.g., Melott 1990) for 
measuring the Euler characteristic (or the genus g = 1 — x) of iso-density contours which 
requires the choice of two parameters, namely a smoothing length and a density level, the 
present method only involves the radius of the covering balls as a single scale. 

• As a consequence of their additivity, the Minkowski functionals are robust and global 
measures: They yield stable results even for small samples. 

• The method is efficient in discriminating different spatial patterns and theoretical models. 
Moreover, the method does not depend on assumptions about the statistics of patterns, or 
model realizations, respectively. 

• The mean values of the Minkowski functionals provide statistically unbiased descriptors 
which contain features of n-point correlation functions at any order n. Therefore, the 
present method is complementary to the conventional approach restricted to lower-order 
correlation functions. 

Finally, the stereological properties of Minkowski functionals (see, e.g., Weil 1983) 
offer means to infer morphological characteristics of three-dimensional samples from lower- 
dimensional sections such as "pencil beams" or two-dimensional all-sky surveys. These 
stereological applications as well as a more detailed statistical analysis is postponed to 
future work. 

Acknowledgements: We would like to thank John A. Peacock for providing us the new 
Abell cluster catalogue, and Jiirgen Ehlers for valuable remarks on the manuscript. TB is sup- 
ported by DFG (Deutsche Forschungsgemeinschaft). 
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Figure Captions 



Figure 1: (a) Parallel body of a convex polygone in d = 2. (b) The cylindrical segment at 
the edge k with dihedral angle a of a convex polyhedra contributes eaL to the parallel volume; 
therefore 3W2(k) = aL. 

Figure 2: Illustrative examples of the Euler characteristic in two and three dimensions. In 
two dimensions \ is the difference of components and holes, whereas in three dimensions it is 
the sum of the disconnected components and cavities minus the number of tunnels. A negative 
Euler characteristic indicates network-like structures. Bottom: Illustrative calculation of \ using 
additivity. 

Figure 3: Mean values of the Minkowski functionals, i.e., of the reduced surface area $i(x) 
(full rectangles), the reduced integral mean curvature $2(3^) (empty rectangles), and the reduced 
Euler characteristic $3 (a;) (triangles) are depicted for a Poisson-process (3a), a 'double-Poisson' 
process (3b) and the ACO-catalogue (3c). The curves in Fig. 3a are obtained from the analytical 
formulas (eqs. (24)). The measures are normalized to 1 at x = and the radius x of the balls 
centered at each point of the process is measured in units of the mean distance p -1 / 3 of the 
points, where p is the number density of the process. 

Figure 4: Thin planar slices from a three-dimensional cube: the upper panel shows 500 points 
which are independently distributed (Poisson-process); the lower panel shows 750 points which 
are distributed according to the 'double-Poisson' process as explained in Section 3.3. 

Figure 5: The mean Euler characteristic per surface area X*{ x ) f° r the three point processes 
discussed: a Poisson-process (P), a 'double-Poisson' process (DP), and the ACO-catalogue (ACO). 
The zeros near x p = 0.44 indicate similar percolation thresholds. The discrepancy of (P) and 
(ACO) obviously rules out the Poisson-process as a statistical model, whereas the similarity of 
(DP) and (ACO) is remarkable. 

Figure 6: Subsamples from the ACO-catalogue with the individual Euler characteristic x °f each 
ball (with reduced radius x) as selection criterion, projected onto a plane. This illustrates the 
possibility to classify cluster environments by using the concept of local morphology introduced 
in this work. 
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